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We present measurements and calculations of the trap loss rate for laser cooled 87 Rb atoms con- 
fined in either a magneto-optic or a magnetic quadrupole trap when exposed to a room temperature 
background gas of Ar. We study the loss rate as a function of trap depth and find that copious 
glancing elastic collisions, which occur in the so-called quantum-diffractive regime and impart very 
little energy to the trapped atoms, result in significant differences in the loss rate for the MOT 
compared to a pure magnetic trap due solely to the difference in potential depth. This finding 
highlights the importance of knowing the trap depth when attempting to infer the total collision 
cross section from measurements of trap loss rates. Moreover, this variation of trap loss rate with 
trap depth can be used to extract information about the differential cross section. 
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I. INTRODUCTION 

It is well known that the lifetime of a trapped ensemble 
of ultracold gas is limited by collisions with atoms and 
molecules present in the residual background vapor of an 
ultra-high vacuum system. The particle loss rates and 
corresponding collision cross sections of laser cooled al- 
kali metals with various species present in a background 
vapor have been studied since the first demonstrations 
of magneto-optic traps (MOTs) [TJ |U [3]. Small angle 
elastic collisions, which impart a very small amount of 
kinetic energy, can result in ensemble heating without 
immediate particle loss. This effect has been observed 
experimentally [4 and analyzed theoretically both in the 
classical and the purely diffractive limits (5j [6l [7j [8] . Elas- 
tic as well as inelastic collisions with a background vapor 
present an array of decoherence channels which can act 
on both the internal and external degrees of freedom of 
a trapped particle. Understanding the exact character 
of these collisions is of fundamental importance for the 
production of ensembles of cold atoms or molecules and 
their use, for example, in quantum information storage 
and processing applications. 

While generally considered a nuisance, particle loss 
from a MOT induced by collisions with a background gas 
was recently proposed as a reliable observable to deter- 
mine the total collision cross section. Matherson et. al re- 
ported measurements of the absolute total collision cross 
section between various room temperature gases (both 
atomic and molecular) and laser-cooled metastable neon 
atoms by detecting their loss rate from a MOT [9j [10] . 
It was argued that trap loss measurements of this sort 
can be done with a very high precision and that they 
are superior to measurements of the cross sections using 
crossed atomic beam experiments where the determina- 
tion of the absolute collision cross section is complicated 



by uncertainties in the absolute number of target atoms 
and the overlap of the beams [TU]. This MOT trap loss 
technique is similar in spirit to previous work where the 
measurement of absolute photoionization cross sections 
of atoms jTTJ [T^] and of absolute electron-impact ion- 
ization O E] were made using laser cooled atoms as 
a target. In the work of Matherson et. al, it is claimed 
that their MOT loss measurements were not significantly 
affected by the retention of collisional products and that 
the cross section for trap loss was therefore equivalent to 
the total collision cross section [10]. By contrast, Steane 
et. al argued that the collision limited trap lifetime of 
a Cs MOT was significantly influenced by the recapture 
of collisional products [3]. Similarly, Gensemer et. al re- 
ported dramatic variations in background collision loss 
rates from a Rb MOT as the trap parameters and depth 
were varied [13]. In the work of Schappe et. al a Rb MOT 
was intentionally operated such that the ions produced 
by electron collisions would escape from the trap while 
atoms that were simply excited or elastically scattered 
were retained [14] . 

In this work we investigate the use of trapped, laser- 
cooled rubidium atoms to measure the cross section for 
87 Rb- 40 Ar collisions. Motivated by the known role of 
small angle scattering in ensemble heating due to col- 
lisions that do not directly result in particle loss, we 
studied the trap loss rate produced by room-temperature 
background collisions as a function of trap depth at 
1 K and below 10 mK using a magneto-optic and a 
quadrupole magnetic trap. We find that for room tem- 
perature Ar, the retention of Rb atoms due to a finite 
trap depth (on the order of or larger than the energy 
transfer associated with quantum diffractive collisions) 
significantly suppresses the trap loss rate in comparison 
with the total collision rate. Using model interaction po- 
tentials for the Rb-Ar complex, we confirm that the ex- 
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perimentally measured loss rates at various trap depths 
are consistent with the known long-range attractive part 
of the potential characterized by the C 6 coefficient. These 
experimental and theoretical results are summarized in 
Fig. [9] and it constitutes the central result of this paper. 
The strong variation of the trap loss rate as a function of 
the trap depth even at depths far below 1 K illustrates 
the very low energy-scale and abundance of quantum- 
diffractive collisions, which account for approximately 
half (exactly half for a hard core potential) of the total 
cross section [B] . The conclusion is that by using a very 
shallow magnetic or optical dipole trap, whose depth can 
be easily characterized and made to approach zero, the 
total loss rate can be made to approach the total collision 
rate. In the context of molecular beam measurements of 
the total collision cross section, a zero trap-depth mea- 
surement is equivalent to realizing an ideal apparatus of 
unlimited angular resolution |16j . 



II. THEORY 

When a room-temperature background vapor particle 
encounters a trapped atom, the collision energy is typi- 
cally very high. Nevertheless, for glancing (very high im- 
pact parameter) elastic collisions, the momentum trans- 
ferred and resulting scattering angle can be extremely 
small. Consequently, for a trap with finite depth, col- 
lisions below a certain scattering angle do not result in 
trap loss. The classical small- angle approximation pre- 
dicts the cross section for loss inducing collisions increases 
without bound as the trap depth is reduced to zero. This 
is because the inter-atomic potential is infinite in range 
and collisions at all ranges of impact parameter impart 
some, albeit small, momentum to the target atom thus 
producing trap loss. However, the total cross section for 
loss inducing collisions can never exceed the total colli- 
sion cross section which is finite and, through the opti- 
cal theorem, related to the imaginary part of the quan- 
tum mechanical scattering amplitude [IT]. In the limit 
of vanishingly small scattering angles, the associated de 
Broglic wavelength of the momentum transfer in the col- 
lision exceeds the classically determined impact parame- 
ter and a classical small-angle approximation is no longer 
valid [18j [19]. In this limit, the collisions are said to 
be quantum-diffractive and a full quantum-mechanical 
treatment of the scattering cross section is needed |20j . 

Consider a trapped atom of velocity v a and mass M a 
and a background particle of velocity and mass M b . 
The initial relative velocity is v r = v a — ilh, and mo- 
mentum conservation requires that the change in these 
velocities is related by — M a Av a = M^Av^ and Av a = 
-^-Av r , where the reduced mass is \i — j^ Ia ^f} b ■ In ad- 
dition, energy conservation requires that for elastic col- 
lisions |Av r | 2 = 2|v r | 2 (l — cos 9), where 9 is the colli- 
sion angle between the initial (v r ) and final (v T ') rel- 
ative velocities. After the collision, the trapped atom 
kinetic energy will have changed by an amount AE = 



[(v a + Av a ) 2 - = ^ [2v a ■ Av a + \Av a \ 2 ] . If this 
change in kinetic energy exceeds the trap depth Uq, loss 
will occur. In the limit that the trapped particle has a 
negligible initial kinetic energy (w a <C Av a ), we have that 
AE~ ^|Av a | 2 and 

AE^^-\v T \ 2 (l-cos9). (1) 

Therefore, if the collision angle exceeds the minimum an- 
gle 

flmin = arccos (l - * ) , (2) 

then trap loss will occur. The rate at which background 
particles are scattered from a single trapped atom into a 
solid angle d£l is n^v^da / d£l)d£l, where is the den- 
sity of background particles and (da/dU) is the differen- 
tial scattering cross section. Given this, we can estimate 
the loss rate from a trap of depth Uq as 

7io SS = ^bWprobfioss = n h v ploh / (da I 'dil,)dn. 

(3) 

where the relative collision velocity is assumed to be de- 
termined by the most probable velocity for the back- 
ground particles, \v r \ = t> pro b- This loss rate is al- 
ways smaller than the estimated total collision rate, 
7c = "-bWprobC, since 9 mm > 0. The differential cross 
section is related to the quantum mechanical scattering 
amplitude, da/dfl = \ f(k, 9)\ 2 , and depends explicitly on 
the collision wave vector, k — {j,\v v \/h. It is assumed that 
the interaction potential is central and therefore there is 
no azimuthal angular dependence. For a beam of inci- 
dent scattering particles with wave vector fc, the cross 
section for loss inducing collisions from a trap of depth 
U is 

a loss (k)= f 2irsm9\f(k,9)\ 2 d9. (4) 

Je min (ftfc/Ai) 

This expression is equivalent to the total cross section 
measured by a molecular beam apparatus with finite an- 
gular resolution limited to (9 m i n [16 . 

To compute the loss rate induced by collisions with a 
background gas at temperature T, we need to average 
over the Maxwell-Boltzman distribution of incident col- 
lision wave vectors. Assuming the trapped atom velocity 
is negligible, we have that v T ~ Vb and we can compute 
the velocity averaged loss rate using k = fi\vb\/h and find 
(Tloss) = n b (va) where 

(v*) - (^) 3/2 /°° 47r ^( fc )-be- Mb ^ /2fcBT ^b. 

(5) 

In order to evaluate this expression, we need to deter- 
mine the scattering amplitude, f(k,9), given the inter- 
atomic potential between atom a and b. The asymptotic 
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form of the 2-body scattering wave function is the sum 
of an incident plane wave and a scattered spherical wave 

"0( r )|r^oo = A (e ikz + f(k, 0)^P^j . Given a central po- 
tential, ip and / can be expanded in terms of the Legendre 
polynomials, 



Mr, 9) = J2Ri(k,r)P l (cost 
1=0 

oo 



(G) 
(7) 



1=0 



This is an expansion into partial waves of angular mo- 
menta I. For sufficiently large r, the potential is negli- 
gible, and the radial functions asymptotically approach 
the form for a free particle 



Rl(k, r) | = Biji(kr) +Cim(kr), 



(8) 



where ji and n\ are the spherical Bessel and Neumann 
functions. The coefficients Bi and C; can be written as 
Bi = Ai cos 5i and Ci = —AisinSi, where Ai = (21 + 
l)i l e lSl and where Si = arctan(— Ci/Bi) is the phase shift 
of the I partial wave [5TJ [55] . These phase shifts also 
determine the scattering amplitude 



oo 

/(M) = r X)(2I + l)e Wl sin ^P, (cost 



(9) 



The determination of the scattering amplitude and re- 
sultant collision cross section therefore requires finding 
the partial wave phase shifts, and they are found by nu- 
merical integration of the radial Schrodinger equation. 
We write the radial equation for the I th partial wave as 



dr 2 



where 



Wi(r) 



2/i 



Mr) = 



Wi{r) = k 2 --^V(r) 
ipi(r) = krRi(r), 



1(1 + 1) 



(10) 



(11) 
(12) 



and V(r) is the inter-atomic potential. The solution to 
the radial equation for each partial wave I is then in- 
dependently computed using the log-derivative method 
[23] . The integration of the radial equation starts from a 
position deep inside the classically forbidden region and 
ends at a point in the asymptotic regime, beyond which 
the contribution becomes insignificant. The phase shift is 
found by matching this numerical solution in the asymp- 
totic region to the asymptotic form, Eq. [8] and by using 

The resultant values for partial wave phase shifts Si de- 
pend on the incident wave vector magnitude k — fi\v r \/h. 
Given the results for a particular collision wave vector, 
the scattering amplitude can be constructed using Eq. [9] 
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FIG. 1: (Color online) Schematic of the vacuum system. The 
variable leak valve allows Ar gas to be introduced and the 
residual gas analyzer (RGA) is used to measure the pressure 
and composition of the background vapor. The laser cooled 
and trapped 87 Rb atoms undergo collisions with the back- 
ground gas which may or may not result in trap loss. 



and inserted into Eq. [4] to compute the total cross sec- 
tion for loss inducing collisions at that collision velocity. 
Repeating this procedure for a set of wave vectors chosen 
from a Maxwell-Boltzman distribution, the velocity av- 
eraged rate of loss inducing collisions from a background 
gas at temperature T (Eq.[5| is evaluated. In each numer- 
ical simulation, the convergence of these results is checked 
both as a function of the integration range and step size 
and as a function of the number of partial waves included. 
It is found empirically that the results (for room temper- 
ature collisions) have converged to an accuracy of better 
than 1% when partial waves up to I — 500 are included 
(higher velocities require more partial waves) and given 
the integration is performed with a step size of 0.005 ae 
(&b = 0.529177 A) from 1 ae, deep inside the classically 
forbidden region, to 150 ae, well in the the asymptotic 
regime. 



III. EXPERIMENTAL SETUP 

Our apparatus consists of a magneto-optic trap (MOT) 
which collects and cools a few million rubidium atoms 
from a room-temperature Rb vapor of below 10 -8 torr. 
The atomic vapor is maintained by periodically heat- 
ing an Alvatec rubidium source. The laser system used 
for 87 Rb laser cooling is composed of grating-stabilized 
and injection-seeded diode lasers and has been described 
previously [21] ■ After amplification, we have a total of 
20 mW of pump or cooling light (driving transitions from 
\F = 2) — > \F' — 3)) and 4 mW of repump light (driv- 
ing transitions from \F = 1) — > \F' = 2)) both oper- 
ating on the 5 2 <Si/2 — ► 5 2 i-3/2 transition. This light 
is combined and then expanded to a 1/e 2 diameter of 
7 mm, prepared with the correct polarization and intro- 
duced into the MOT cell (a borosilicate glass box) along 
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the three mutually orthogonal axes in a retro-reflection 
configuration. We typically operate with an axial gra- 
dient of b' — 23.1(0.7) Gem" 1 and with pump and re- 
pump dctunings of Ai/ p = — 2"/m> and Az/ r = 0, where 
7Rb = 2-7T x 6.07 MHz is the natural linewidth of the 
5 2 <Si/2 - * 5 2 P3/2 transition. The total pump intensity is 
I = 5.6/ sat (where I sat = 3.58mWcm~ 2 for an isotropic 
pump field with equal components in all three possible 
polarizations |25jl. Based on previous measurements of 
the escape probability as a function of energy for a Rb 
MOT [JB] and studies of the dependence of the capture 
velocity on intensity jTB] HZ] , we estimate the MOT trap 
depth in this work to be 800 ± 300 mK. The number of 
trapped atoms in the MOT is determined by measuring 
the emitted resonance fluorescence. 
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FIG. 2: (Color online) MOT Fluorescence as a function of 
time fit to the rate equation given in Eq. [13] with f3 = 0. 

87 Rb atoms are cooled, optically pumped into cither 
the \F = 1) or \F = 2) hypcrfine ground state and 
then transferred from the MOT into a quadrupole mag- 
netic trap. For this purpose, the MOT anti-hclmholtz 
coil pair is used and provides an axial gradient of up 
to b' — 670(20) Gem -1 . The intrinsic loss rate in 
this quadrupole magnetic trap due to Majorana spin 
flips is negligible at the temperatures and densities in 
this experiment |28j . The atom number remaining in 
the magnetic trap after some hold time is measured by 
rapidly returning the magnetic field to an axial gradient 
of b' = 23.1(0.7) Gem -1 , used for the MOT, switching 
the MOT lasers on and detecting the fluorescence of the 
atoms which are recaptured in the MOT. Within a mil- 
lisecond, the atoms are laser cooled and recollected in 
the MOT and the resonance fluorescence is monitored 
during the first 100 ms after recapture. The measured 
atom number is determined from this fluorescence and 
is corrected for the number of atoms captured from the 
residual 87 Rb background vapor during the measurement 
time. Since the trap loss rates reported here are deter- 
mined from the ratios of atom numbers remaining in the 
trap after various hold times, the only requirement for 
a precise determination of the loss rate is that the mea- 
sured fluorescence be linear in the number of trapped 
atoms. Since the optical density is well below 0.5 for all 
of the measurements reported here, this is the case. 

Additional room-temperature Ar gas is introduced into 
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FIG. 3: (Color online) Fraction of atoms retained in a 
quadrupole magnetic trap as a function of time fit to the rate 
equation given in Eq. [13] with R = and /3 = 0. 

the vacuum envelope through a variable leak valve. This 
gas is continuously pumped away by a turbo-molecular 
pump, and the resulting equilibrium partial pressure is 
controlled by the leak rate of the value and measured 
using a residual gas analyzer (SRS RGA 200) as shown 
in Fig. [T] The additional trap loss rate induced by this 
new gas as a function of its pressure is measured. 

Although an accurate determination of the exact cross 
section using this trap loss method would require care- 
ful calibration of the residual gas pressure, the focus of 
this work is to illustrate the dependence of the collision 
cross section on trap depth. Therefore this work requires 
only that the the ratio of cross sections at different trap 
depths be precise. In order to ensure this, and to mitigate 
the effect of the uncertainty in the absolute pressure, the 
data were collected as follows: the Ar pressure was first 
held fixed and the loss rates for the MOT and magnetic 
trap at various depths were measured. The pressure was 
then adjusted and the procedure was repeated. This en- 
sured that although the absolute pressure reading may 
be in error, the loss rates in traps of varying trap depth 
were all measured at the same pressures. The ratio of 
cross sections between traps of various depth is therefore 
precise to within a statistical error of less than 10%. 

In the MOT, the atom number loading dynamics in 
the presence of a background vapour is modelled by the 
following rate equation [2"5] , 

^=R-TN-pJn 2 Rh dV (13) 

where N is the trapped atom number, R is the load- 
ing rate, T is the loss rate due to collisions between 
the trapped rubidium atoms and the room-temperature 
background gas, and (3 is the combined density depen- 
dent 2-body loss coefficient due to cold collisions within 
the MOT. Here 3-body and higher order collisions are 
neglected. The density profile of the cold, trapped atoms 
is given by n^b, and these density dependent losses in- 
clude light-assisted collisional losses (such as radiative 
escape and fine structure changing collisions) and hypcr- 
fine structure changing collisions. For sufficiently small 
atom numbers, the last term is negligible compared to the 
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second [24 . In the case of the pure magnetic quadrupole 
trap, the trapped number is well described by this model 
with R — 0. In the magnetic trap, the confinement is 
much weaker than the MOT and the trap volume (up 
to 5 mm in radius compared to the MOT cloud radius 
of typically 150 fim) is larger by a factor of more than 
10 4 , making the density dependent collisional losses even 
smaller. In addition, light assisted collisions (which dom- 
inate the j3 term for the MOT) are no longer present, and 
only hyperfine changing collisions and spin relaxation col- 
lisions (from magnetic dipolc-dipole coupling) remain. In 
any case, the data presented here are taken in the regime 
where these 2-body loss processes are completely domi- 
nated by the TN term. 



IV. RESULTS 

Figures[2]and[3]show the loading dynamics of the MOT 
and the atom loss from a magnetic trap fit to Eq. [13] 
From these fits, T is extracted. The loss rate due to 
collisions with the room temperature gas includes a con- 
tribution from three different terms 

r = r^b + r rcsi duai + Tb ■ (14) 

r^b is the loss rate due to collisions with room- 
temperature Rb atoms, and r rcsic i U ai is the loss rate due 
to collisions with unwanted residual gas in the ultra-high 
vacuum. Tb is the loss rate due to collisions with gas 
particles intentionally introduced into the vacuum sys- 
tem and is given by T h = n h (a Rhyh v h ) where (cr Rbib Wb) 
is the velocity averaged collision cross section of interest 
(Eq.|5f. 

As the controlled background gas density is in- 
creased, the total loss rate, T, increases with slope 
(cRb.bWb)- In general, all three of the loss rates con- 
tributing to r (r^b, Tb and r rcs j(3uai) can be different in 
the MOT than in a magnetic trap. This is due to the 
difference in trap depths and to the fact that some of the 
atoms in the MOT are in an electronic excited state dur- 
ing a collision (on the order of 12% for our MOT param- 
eters) which will, in general, have a different interaction 
potential and corresponding cross section with the collid- 
ing particles than do ground state Rb atoms. The case 
of excited state Rb atoms colliding with identical ground 
state Rb atoms is an extreme example. The leading order 
in Rb-Rb ground state interactions is Ce/r 6 whereas it is 
C^/r 3 for Rb-Rb* collisions. In this way, a measurement 
of magnetic trap loss is less complicated to interpret since 
it only involves elastic collisions of ground state atoms. 

Figure |4] shows the measured loss rate of trapped 87 Rb 
versus argon density in both a MOT and a quadrupole 
magnetic trap. These data are well fit by a line and 
the slopes, which provide the value of the velocity aver- 
aged collision rate (<7Rb,Ar ^Ar), are listed in table III The 
MOT exhibits the smallest cross section (i.e. slope) due 
to a much larger trap depth (~1 K) than the magnetic 
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FIG. 4: (Color online) Loss rate of trapped 87 Rb versus room 
temperature argon gas density. The results are fit to a line and 
the slope provides the value of the velocity averaged collision 
rate (o"Rb,Ar Vax). The MOT (circles) exhibits the smallest 
cross section due to a much larger effective trap depth (~1 K) 
than the magnetic trap (squares) where the \F — 1, mF = — 1) 
state was confined with b' = 79 Gem" 1 resulting in a trap 
depth of 0.64 mK. The vertical error bars, most of which are 
smaller than the symbols, represent the statistical uncertainty 
in the loss rate based on fits similar to those shown in Figs. [2] 
and [3] 



trap (~1 mK). As the argon density and associated back- 
ground collision loss rate are increased, the steady state 
MOT number is smaller as are any density dependent 
losses (our MOT is operating in the small-particle range 
where the density increases linearly with atom number 
[2~4"]y If the density dependent losses were significant, 
one would therefore expect the loss rate to vary non- 
linearly with argon gas density. The MOT loss rate is 
observed to vary linearly with the argon density confirm- 
ing that the combined density dependent 2-body losses 
due to cold collisions within the MOT are negligible un- 
der our experimental conditions. 

Before loading the pure magnetic trap, hyperfine opti- 
cal pumping is employed to populate and confine either 
the \F = 2) or the \F = 1) state. For small magnetic 
fields where the electronic and nuclear magnetic moments 
are strongly coupled by the hyperfine interactions, F is 
a good quantum number and the \F — I, nip = — 1) 
and \F = 2, tuf — 1, 2) states arc the magnetically trap- 
pable, weak-field seeking states. Also at small magnetic 
fields, the |F = l,m F = -1) and \F = 2,m F = 1) 
states have essentially identical magnetic moments (dif- 
fering by less than 0.0004% at zero magnetic field) while 
the \F = 2,rriF = 2) state has twice the magnetic mo- 
ment of the others and therefore experiences a larger 
trap depth for a given magnetic field gradient. The 
quadrupole coils are oriented with the axis of symme- 
try along the direction of gravity (with an axial gradient 
of up to b' = 670(20) Gem -1 and a radial gradient of 
b' /2). For small fields, the magnetic trap potential en- 
ergy along the axial direction (z) is £/g xlal = ab'\z\ + mgz, 
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TABLE I: Measurements of the trap loss rate slope for 87 Rb 
in the presence of room temperature 40 Ar gas listed in order 
of increasing trap depth. These results are compared with 
the theoretical computed trap loss rate slope shown in Fig. [9] 
The minimum and maximum values within the uncertainty in 
the magnetic trap depth indicate the depth along the radial 
and axial directions while the uncertainty in the MOT trap 
depth represents the estimated trap depth range based on the 
work of Hoffmann el. al 1261. 



Trapped State 
\F,m F ) 


b' 

(Gcm^ 1 ) 


Trap Depth 
(mK) 


(<TR,b,Ar VAt) 

(xl0- 9 cm 3 s _1 ) 


|1,-1> 


39.5 


0.24 (0.09) 


2.45 (0.07) 


|2,+2) 


23.7 


0.34 (0.06) 


2.40 (0.11) 


|1,-1> 


79.0 


0.73 (0.07) 


2.20 (0.09) 


|1,-1> 


111 


1.13 (0.20) 


2.18 (0.12) 


|1,-1> 


334 


3.77 (1.04) 


1.98 (0.09) 


|1,-1> 


501 


5.64 (1.59) 


1.86 (0.06) 


|1,-1> 


668 


7.41 (2.09) 


1.64 (0.04) 


MOT 


23.7 


800 (300) 


0.54 (0.02) 



where a — fiBgF'm-F is the magnetic moment of the 
state, fiB — h x 1.4 MHz/G is the Bohr magneton, and 
g F = l/2(— 1/2) for the \F = 2) (\F = 1)) state. Along 
the radial direction, the trap potential energy is, for small 
magnetic fields, [/o' adlal = a\r. 

For large magnetic fields, the electronic and nuclear 
spin magnetic moments begin to become decoupled and 
the Zecman shift is non- linear. The trap depths reported 
(Tableland Fig.|J were evaluated by diagonalizing the 
full Hamiltonian including the hyperfine and Zeeman in- 
teractions. The magnetic trap depth is limited by colli- 
sions with the cell walls which are 5 mm from the mag- 
netic trap center along the radial and axial directions. 
Since particles can escape along either the axial or radial 
directions, a trap depth error bar is provided in Table [T] 
and Fig. [9] indicating the depths along these directions. 
Typically, the trap depth is limited by radial trajecto- 
ries; however, at axial gradients below b' < 2 ^ gravity 
limits the trap depth along the z direction. In addition, 
below a critical gradient of b' < ^ = 30.5 Gem -1 , only 
the \F = 2, mp = 2) is supported against gravity. This 
feature can be used to "clean" the \F — 2) manifold re- 
moving all atoms in the other low-field seeking states. As 
discussed above, for the lower hyperfine state, only the 
\F = X,mp = —1) state is magnetically trapped, while 
for the upper hyperfine state, due to the non-linear Zee- 
man shift, there are three states (which correlate to the 
\F = 2, mp — 2, 1, 0) states at zero magnetic field) which 
can be trapped magnetically, and each state has a dif- 
ferent trap depth. To simplify the interpretation of the 
data, we focus on the loss rate of the \F — l,mp = —1) 
state. In addition, we study trap loss for the upper hy- 
perfine state at a field gradient of 23.7 Gem -1 . At this 
gradient only the \F = 2,rriF = 2) state is trapped and 
this ensures a single trap depth for the trapped atomic 
ensemble. 

In this analysis, the inter-atomic interaction is modeled 
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FIG. 5: (Color online) The theoretically computed total cross 
section for room temperature 87 Rb- 40 Ar collisions (averaged 
over a Maxwell-Boltzmann distribution at 295 K) given a 
Lennard-Jones inter-atomic potential as a function of the Ce 
coefficient (C12 = 8.6 x 10 7 E^a^). The numerical results are 
fit to a power law, a — 70A 2 x (Ce) ' °, in agreement with 
the approximate dependence of 2/5. 



by a Lennard-Jones potential, V(r) = Cyijr 12 — Cq/v 6 , 
where the value Cq = 280 E^a^ is used for the Rb-Ar 
interaction (E h = 4.35974 x 10" 18 J). This coefficient is 
known from both optical and atomic beam measurements 
[501 EH E2- C12 = 8.6 x 10 7 E h a£ is chosen so that 
the resulting potential has a minimum of 50 cm -1 below 
zero. This is the order of magnitude expectation for the 
Ar-Rb complex [33]. However, this choice is somewhat 
arbitrary, since for the collision energies involved in this 
experiment, the exact value of C12 has a negligible effect 
on the total cross section. For a long-range potential of 
the form V(r) = —C n /r n , where n > 2, it can be shown 
[T7] that the cross section can be expressed as 



o-(fc) = A 



9n 
hv 



2/(n-l) 



(15) 



where A is a constant and v is the collision velocity. This 
expression is approximate and only valid for sufficiently 
large values of the collision wave vector /iC n k n ~ 2 / ti 2 ^> 1 
(satisfied for all Ar-Rb collision velocities above 1 m/s 
for Cq — 280 E^a B ) and for potentials which have the 
form specified for distances beyond r c ~ {C n /hv) l /^ n ~ x ' . 
Given the known value for C§ and our choice for C12, 
the Lennard-Jones potential deviates from a pure Ce /r 6 
potential by less than 3% at distances beyond r c , defined 
for Ar-Rb collision velocities of 1000 m/s. The devia- 
tion is even smaller for lower collision velocities since r c 
is larger. The result is that for collisions with a room 
temperature gas, the cross section is insensitive to the 
exact value of C12. Of course, for much higher velocities, 
the exact short range behavior of the inter-atomic poten- 
tial begins to matter, and at velocities much lower than 
1 m/s, low partial wave collisions are dominant and the 
above approximate expression is invalid. 

As a check of the numerical calculations, the computed 
cross section is compared with the predicted power law 
dependence on Cq, and an exponent of 0.40 is found, 
in agreement with the expected 2/5 dependence from 
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Eq. 15 The results are shown in Fig. [5] 
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FIG. 6: (Color online) The theoretically computed total cross 
section for 87 Rb- 40 Ar collisions averaged over a Maxwell- 
Boltzmann distribution at various temperatures given a 
Lennard- Jones inter-atomic potential with Cg = 280 E^a^ 
and C12 = 8.6 x 10 7 E^ajj. The results are fit to a power 
law, a = 2100A 2 x (j 1 )-° 202 i m agreement with the approx- 
imate expected power dependence of -1/5, corresponding to 
the power dependence of -2/5 on the collision velocity. 

The total cross section for 87 Rb- 40 Ar collisions aver- 
aged over a Maxwell-Boltzmann distribution, (<7Rb,Ar), 
is also computed at various Ar temperatures and shown 
in Fig.[6l These results also agree with the predictions 
of Eq. |15[ In addition, they are consistent with previ- 
ous atomic beam measurements of Ar-Rb collisions. In 
that work, the attenuation of a 283 K Ar beam travel- 
ing through a 590 K Rb vapor was measured and a cross 
section of 572 A 2 was found |30J. Although the average 
beam velocity was not specified, we estimate that the 
relative collision speed in that measurement corresponds 
to a stationary Rb atom exposed to a T > 400 K Ar 
gas. Our calculations predict a cross section of 626 A 2 at 
T = 400 K, within 10% of the value reported in Ref. [3D]. 
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FIG. 7: (Color online) The theoretically computed total 
cross section for room temperature 87 Rb- 40 Ar collisions (av- 
eraged over a Maxwell-Boltzmann distribution at 295 K) 
given a Lennard- Jones inter-atomic potential as a function 
of the C12 coefficient with Ce = 280 E^a%. A value of 
Cia = 8.6 x 10 7 £ h a B 2 results in a potential with a minimum 
of 50 cm -1 below threshold. As expected, for the collision 
energies characteristic of a room temperature gas, the exact 
value of the G\i coefficient (at or below 10 9 E^ajj) plays little 
role in the resulting cross section. 

Finally, the sensitivity of the exact value of the C\% co- 



efficient on the computed cross section was also checked. 
The results of those calculations are shown in Fig. [7] The 
conclusion is clear, for any values of C12 less than or equal 
to 10 10 Eh^e 2 ' more than 10 times larger than our cho- 
sen value, the cross section remains essentially unchanged 
(varying by less than 3%). 

In Fig. [8] the theoretically computed total loss cross 
section of trapped 87 Rb as a function of the trap depth 
for 87 Rb- 40 Ar collisions at different collision velocities is 
shown. As one would expect, the variation in the loss 
inducing cross section with trap depth is the most pro- 
nounced at low collision velocities. However, there still 
remains a significant variation in cross section for a room 
temperature thermal distribution of collision velocities. 




10 
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FIG. 8: (Color online) The theoretically computed loss cross 
section as a function of trap depth for 87 Rb- 40 Ar collisions 
of different velocities. Loss cross sections for velocities corre- 
sponding to the average speed at 4 K (open triangles), 295 K 
(open diamonds), 1000 K (open squares) are shown as well as 
the average loss cross section for a room temperature thermal 
distribution (filled circles). The variation of the loss cross 
section with trap depth is most pronounced for the lowest 
velocity collisions. However, there still remains significant 
variation with trap depth for a thermal distribution at room 
temperature. 

Fig. [9] shows the theoretically computed and exper- 
imentally measured loss rate slope, (oRb.Ar i>Ar), for 
trapped 87 Rb as a function of the trap depth. This is 
the central result of this paper. It reveals that the trap 
loss rate is a strongly varying function of the depth down 
to values below 1 mK. In fact, the variation with depth 
is observed to become more pronounced at values just 
below 10 mK. This behavior can be understood in the 
following way: approximately half of the total collision 
cross section arises from classical scattering (exactly half 
for a hard core potential) while the other half arises from 
quantum diffractive collisions |21j . The natural energy 
scale for diffractive collisions is 



AttH 2 



(16) 



where M a is the trapped atom mass and a is the to- 
tal collision cross section [6j. For Ar-Rb collisions, this 
energy scale is 10 mK, and a trap of this depth would 
exhibit roughly half the loss rate of a trap of zero depth. 
The computed loss rate slope at 10 mK is (<7Rb,Ar VAi-) — 
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1.27 x 10~ 9 cm 3 s _1 , almost exactly half the zero trap 
depth loss rate slope of 2.43 x 10 -9 cm 3 s _1 . For trap 
depths below 0.01 e^, the Rb loss rate differs from the 
zero trap depth rate by less than 1%. The conclusion is 
that by using a very shallow trap, whose depth can be 
easily characterized and made small compared to ea, the 
total loss rate can be made to approach the total collision 
rate. We note that the variation of the loss rates within 
the range of experimentally accessible trap depths pro- 
vides a means to map out the differential cross section 
from the diffractive regime to the classical small angle 
regime. Measurements of this sort will be the subject of 
future work. 
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V. CONCLUSIONS 



In summary, the use of trapped, laser cooled rubidium 
atoms to measure the cross section for 87 Rb- 40 Ar colli- 
sions is investigated. Trap loss rates produced by back- 
ground collisions as a function of trap depth at 1 K and 
below 10 mK using a magneto-optic and a quadrupole 
magnetic trap are studied. The retention of atoms due to 
a finite trap depth (i.e. larger than the quantum diffrac- 
tive collision energy scale) is found to significantly reduce 
the measured loss rate in comparison with the total col- 
lision rate. Using model interaction potentials, the ex- 
perimentally measured loss rates at various trap depths 
are found to be consistent with the known long range 
Cq coefficient. These results highlight the importance of 
minimizing the trap depth when attempting to infer the 
total collision cross section from measurements of atom 
trap loss rates. Finally, this analysis could be used to 
accurately compute heating rates in magnetic traps pro- 
viding a complementary approach to analytic methods 

MM HIE]- 



FIG. 9: (Color online) The experimentally measured 
(squares) and theoretically computed (circles) loss rate slope, 
(o"Rb,Ar vai), as a function of trap depth for 8r Rb- 40 Ar colli- 
sions. The line is only a guide to the eye. The data point at 
1 K was obtained using a MOT and the data below 10 mK 
were obtained with a quadrupole magnetic trap. These values 
were extracted from the data shown in Fig [4] The horizontal 
error bar for the MOT data point represents an estimated trap 
depth range based on the work of Hoffmann et. al . The 
error bars for the data below 10 mK represent the radial and 
axial magnetic trap depths. The horizontal dotted line indi- 
cates the zero trap depth limit for (o"Rb,Ar vat), corresponding 
to the total collision cross section. 
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